Predicted light scattering by Josephson vortices in layered cuprates 
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Light scattering by Josephson vortices is theoretically predicted. An experimental set-up for the 
observation of the effect is proposed. The magnetic field should be parallel to the CuC>2 layers and 
the crystal surface. The falling light should be at grazing angles, and the Bragg scattering must 
be almost backscattering. The distance between the columns of the Josephson vortices must be 
slightly smaller than the light-wavelength. The surface should be very clean and flat in order light 
scattering by surface terraces and surface roughness to give the relatively small background. 

PACS numbers: 71.23.An, 72.15.Rn, 72.25.Ba, 73.63.Nm 



I. INTRODUCTION 

The purpose of the present work is to attract the atten- 
tion of the experimentalists to the possible existence of a 
new effect in the physics of the layered high-T c cuprates 
- the light scattering by the Josephson lattice as by a 
grating made from a non-transparent material. 

II. MODEL 

For the description of a typical layered superconduc- 
tor like the high-temperature superconductor (HTSC) 
YBa2Cu307 we use the Lawrence-Doniach (LD) modet^ 
of the Josephson-couplcd superconducting layers. In this 
model the material is treated as a stack of superconduct- 
ing planes, separated by an insulating material. When 
the spacing between the layers tends to zero, the effects 
of the microscopic layered structure are averaged out, 
and the LD model reduces to the anisotropic Ginzburg- 
Landau model. In the latter the anisotropic nature of the 
material appears only in the form of a mass tensor with 
unequal principal values. 

The LD Gibbs free energy of the superconducting lay- 
ers parallel to the xy plane and stacked along the z-axis 
is given by 

^ LD = co / d 2 p F ™ 
n=0,±l,.. 

+ / — [rot A(r) - B(r)} 2 dxdydz, 

where cq is the spacing between the superconducting lay- 
ers, p = (x,y) is the two-dimensional in-plane radius- 
vector. The superconducting layers are enumerated by an 
integer n, with the coordinate of the n-th layer z n = ucq 
(see, for instance, Fig. 4 in Part II of Rcf. 0]) 



arc the three-dimensional vector-potential and space- 
vector. B(r) is the external magnetic field. 

In the representation when the order parameter is the 
wave function of the superconducting electrons, the free 
energy density F n of the n-th superconducting layer from 
the LD functional Eq. (p} is given by 



F„ = a(T)|* n (p)| 2 + ffl|* n (p)| 4 
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where cv(T) and /3(T) are Ginzburg-Landau parameters, 
which are constant in space but depend on temperature, 
^> n (p) and 9 n (p) are the wave function and the supercon- 
ducting phase only defined on the n-th plane, M and m* 
are the effective masses of the Cooper pairs for the motion 
along the z-direction and in the xy-plane (the coordinate 
axes were chosen in such a way, that they coincide with 
the principal directions of the mass tensor of the Cooper 
pairs), e* = 2e is the charge of the Cooper pairs, V 2D 
is the 2D gradient in the xy plane, A n (p) = A{p, nco), 
n = 0, ±1, .., is the restriction of A to the n-th supercon- 
ducting plane. The presence of the vector-potential in 
the free energy, Eq. ([3]), is due to the gauge invariance, 
the phase of the order parameter has to come in the com- 
bination 8 — e* /he J A ■ dr. Thus enclosed in the round 
brackets in the third line of Eq. (|3|) is the gauge-invariant 
phase difference. This is also the discrete analogue of the 
corresponding gradient term. 

The material is isotropic and homogeneous in planes 
parallel to the superconducting layers. However, in the 
perpendicular direction the superconducting behavior is 
governed by the anisotropy parameter 
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The third line in Eq. (|3]) describes the Joscphson in- 
teraction between the layers. We consider the case of 
a extremely anisotropic superconductor £ c <C Co, where 
£,c(T) = fi/ \j2m*\a{T)\ is the coherence length trans- 
verse to the layers direction (when it exists in the con- 
tinuous limit of Co — > 0). In this case, since the center 
of a vortex parallel to the layers is located in the mid- 
dle of the two adjacent planes, both SCsuperconduct- 
ing planes are far from the normal vortex core of size 
£ c in the transverse direction, which is supposed to exist 
in the continuous limit. Therefore, in the third line of 
the Josephson interlayer coupling term Eq. ([3]), for the 
Cooper pair wave function was substituted its equilib- 
rium value ^n(p) = * GL e l6 '' 1 ^, where the modulus is 




(5) 



and roughly equals the average superconducting density. 

Variation of the LD functional, Eq. (TTJ), with respect to 
the parallel and the transverse components of the vector- 
potential 8F/5A n — and HT j&A z = , gives the in- 
planc current 
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with the total current in the perpendicular to the in- 
direction given by j xy = cq 

J2n S ( z ~ nc o)jn, and the 
Josephson current between the n-th and the (n + l)-st 
layers 
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where jj = (c/ 'Ai:\ 2 ){hc/ e*) / 'c is the Josephson critical 
current, and X a b = y m*c 2 /4irne* 2 and A c = X a bT are 
the London penetration depths {n w f GI ). 

Variation of the LD functional, Eq. (TTJ), with respect 
to the phase 8F/56 n = 0, gives the charge conservation 

law V 2D j n + {j^ n+1 - J'" -1 '") /co = 0. Using the latter 

and Eqs. ([6j[7]), one can obtain the vortex-conditions for 
the gauge invariant phase difference. A typical Josephson 
vortex in a layered compound is shown in Fig. 19 in Part 
II of Ref. Q. The magnetic field is along the y-axis. The 
vortex is flattened out along the anisotropy z-axis. The 
structure of the vortex depends on the relation between 
the London penetration depths A a b,A c and the so-called 
Josephson screening length Aj = cqF. 



and correspondingly the Lawrencc-Doniach functional, 
Eq. ([1}, reduces to the anisotropic Ginzburg-Landau, or 
effective mass, functional. Variation of the anisotropic 
Ginzburg-Landau functional (at a zero external field and 
a zero vector-potential) gives the static Ginzburg-Landau 
equation for uniaxial superconductor 
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is the continuous limit of the discrete second derivative. 



Isotropization for extremely type-II 
superconductors 



After the rescaling 
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with the anisotropy parameter F from Eq. ((4]), the 
anisotropic Eq. ([8) transforms to the isotropic Ginzburg- 
Landau equation 
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with the Laplacian 
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A = d 2 x + d 2 + dl 



and \& OL being the modulus of the superconducting order 
parameter for vortexfree space homogeneous phase by 
Eq. ([5]) and £ c is the coherence length in the z-direction. 

The solution of Eq. (jTU)) in the case of a vortex filament 
with an axis parallel to the y-axis is well known 



\J x 2 + z 2 
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(11) 



where Vx 2 + z 2 is the distance from the axis of the fila- 
ment, cj) = arctan(S/a;) is the polar angle round the axis, 
and f(rf) is found from the equation 



A. Continuous approximation 

In the continuous limit of cq — > the expression in the 
square brackets in the third line of Eq. §3§ transforms to 
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The plot of the function f(rf) is shown in Fig. 4 in para- 
graph 30 of Ref. /(0) = and /(oo) = 1. Thus, in 
the continuous limit and in the curved space (x,z), the 
vortex is round and the coherence length £ c is the radius 
of its normal core. 
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However, in our space (x, z) the vortex is oblate. Also, 
one has to take into account the discrete nature of the 
material. For a qualitative evaluation of the effect consid- 
ered in this Letter, we use a trial function which is exact 
for the isotropic extremely type-II superconductors. To 
obtain this trial function, one substitutes the rescaled 
discrete coordinates z = (n — l/2)cor of the planes in 
the solution pip of the Ginzburg-Landau equation for 
isotropic superconductor 



v n (x) = R n (\x\y e ^ x \ 



(13) 



where R n (\x\) = * OL /( y/x 2 + (n - l/^f 5 /^) with 
f{rf) defined by Eq. (JT2J), and 
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being the non-homogeneous superconducting phase in 
the n-th plane. 



III. ELASTIC LIGHT SCATTERING 

Now let us analyze the Josephson current induced by 
the high-frequency vector-potential of the scattered light. 
The static vector-potential of the Josephson vortex has 
a negligible influence, and the Josephson formula for the 
current Eq. ([7]) gives 



where is the z-component of the vector-potential of 
the light 
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polarized in a plane perpendicular to the axis of the vor- 
tex filament 
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For a vortex-free phase 9 n = const we arrive at the 
London current response jz = — (c/4:iT\ 2 )Az n ^ ■ Sub- 
tracting this homogeneous response, we derive the net 
contribution from the Josephson vortex 
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The substitution here of the trial function for the phase 
distribution, Eq. (Ti"4]) . after some algebra gives 

Si^fx)^— 2(c r/2) 2 (in) 

° Jz [ X) 4ir\ 2 c (c T/2) 2 + x 2 Az (U) 

for the current between layers n = and n = 1. The vor- 
tex "core" is "sandwiched" between those layers. Sup- 
posing that the vector-potential is approximately homo- 
geneous at distances of the order of cqT, the integration 



of Eq. flTSJ) with respect to x (using -^rf^ = § ) gives 
the 5-like approximation 
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where AS = itcqT. This formula has a very simple in- 
terpretation a single vortex operates as a small tunnel 
parallel to the y-direction of the external magnetic field, 
with a cross-section AS. 

Now, summing the influence of all the vortices from 
the Josephson lattice 
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This lattice is like an extended Abrikosov lattice along 
the x-dircction. For a strong anisotropy A c /A a 6 = T ^> 1 
the extended lattice of triangles looks like vortex columns 
spaced by distance 
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where $o = ^irhc/e* is the flux quantum. In vertical di- 
rection the distance between two adjacent vortices in one 
column is b — \/2>a/T 2 . This distance is much smaller 
than the wavelength, and if we wish to analyse the light 
scattering, we have to average the current response over 
the vortices in one column. A single vortex operates as 
a hole with an area AS, and a vortex column is equiva- 
lent to a slit with width w = AS/b. In such a way our 
electrodynamic problem is to calculate the intensity of 
light diffracted by regular grating slits with width w. In 
its depths the material is not transparent, that is why 
the light can see only the top end of the vortex columns 
with a height H < Aii gn t- For an order of magnitude 
evaluation of the effect one can take H ~ 100 nm. Along 
this effective depth of the slits one can have many vor- 
tices. At high enough magnetic field the distance between 
the Josephson vortices can be of the order of the lattice 
constant Co, i.e. one vortex between every double CuC>2 
plane. Finally, the top of the columns works as a grating 
with an area A sa S-AS. 

b 

Diffraction of light from such a Josephson vortex lat- 
tice will appear effectively as a diffraction from a grating 
made from a non-transparent material. Our electrody- 
namic problem is to calculate the vector-potential cre- 
ated by the Josephson current, Eq. ([20]) . To extract the 
contribution in the sum responsible for the creation of 
the n-th diffraction maximum in the reflected light we 
expand Eq. (|20[) in Fourier series 



lz{x) = /^CXp <^ 
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where / = (AS/AttX 2 ) (H/b)(c/a), because we substi- 
tuted AS(H/b) for AS in Eq. fl2DJ) (see the reasoning 
in the previous paragraph). 
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The vector-potential created by the n-th current term 
in the Fourier sum Eq. (|22p . which corresponds to the 
n-th diffraction maximum, is 

/t(out) _ 27r£ iqx+iK z—iut /i( in ) (00\ 



where K = ^ - q 2 and q = k { j n) + ^-n. 

From Eq. (|22|) . and considering the light penetra- 
tion depth to be of the order of the light wavelength, 
H ~ Aught, one finds the reflection coefficient for the ra- 
th diffraction maximum, consisting of three factors 




The first factor in Eq. (f2"3")l is intrinsic for the material. 
For YBa2Cu307 with a spacing between CuC>2 planes 
Co ~ 10 A, an anisotropy parameter T sa 30 and a London 
penetration depth in z-direction A c sa 1800 A, it is 




The second factor in Eq. (|23j) is the ratio of the light 
wavelength and the Josephson vortex lattice constant, 
and is adjusted by the external parameters the frequency 
of the falling light and the applied magnetic field. It 
should be of the order of unity, diffraction is favorable 
when the wavelength of the scattered light is compat- 
ible with the spacing between the columns of vortices 
Aught ~ a > ttcqT. The third factor / (ra,6K in ), 0(° ut )) 



is a function of the incident and the reflected angle of 
the n-th diffraction maximum, and the effect can be best 
observed when it is of the order of unity, when the angle 
0( m ' of the falling light is grazing, and 0( out ) corresponds 
to almost backscattering. 

Actually the microscopic calculation using the Green's 
functions methods (sec for example the work by 
Pokrovsky and Pokrovsky^ and the exercise book by Lev- 
itov and ShytovS) gives a (A/hoj) 2 correction and loga- 
rithmic factors which are practically constants (here u is 
the frequency of the falling light and A is the supercon- 
ducting gap). This decreases the magnitude of the ef- 
fect, but we still have enough sensitivity to see the broad 
diffraction maximum. 



IV. CONCLUSIONS 

For a qualitative consideration of the effect we used 
the static approximation. A derivation by means of the 
anomalous Green's functions reduces the diffraction re- 
flection amplitude by (A/hui) 2 , but nevertheless the ef- 
fect survives, and the intensity is not difficult to mea- 
sure. This result should encourage experimentalists to 
consider the perspectives of such a simple experiment for 
the observation of a new phenomenon - light scattering 
by Josephson vortices. 
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